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Abstract. We introduce a notion of resultant of two meromorphic 
functions on a compact Riemann surface and demonstrate its usefulness 
in several respects. For example, we exhibit several integral formulas 
for the resultant, relate it to potential theory and give explicit formulas 
for the algebraic dependence between two meromorphic functions on a 
compact Riemann surface. As a particular application, the exponential 
^ (— i transform of a quadrature domain in the complex plane is expressed in 

terms of the resultant of two meromorphic functions on the Schottky 
double of the domain. 



t> ' 1. Introduction 

. A bounded domain Q in the complex plane is called a (classical) quad- 

rature domain [1], [42], [46], [26] or, in a different terminology, an algebraic 
domain [53] , if there exist finitely many points zi € and coefficients Cf.j € C 
(i = 1, . . . , N, say) such that 



O; hdxdy = Y,Y. c ^ h{3 ' 1] ^k) C 1 ) 

k=l j=l 

^ ■ for every integrable analytic function h in Q [39]. In the last two decades 

there has been a growing interest in the applications of quadrature domains 
to various problems in mathematics and theoretical physics, ranging from 
Laplacian growth to integrable systems and string theory (see recent articles 
[25]), [30], and the references therein). 

One of the most intriguing properties of quadrature domains is their 
algebraicity [1], [22]: the boundary of a quadrature domain is (modulo 
finitely many points) the full real section of an algebraic curve: 

dQ = {z € C : Q(z,z) = 0}, (2) 

where Q(z,w) is an irreducible Hermitian polynomial. Moreover, the corre- 
sponding full algebraic curve (essentially {(z,w) £ C 2 : Q(z,w) = 0}) can 
be naturally identified with the Schottky double f2 of 17 by means of the 
Schwarz function S(z) of d£l. The latter satisfies S(z) = z on d£l and is, in 
the case of a quadrature domain, meromorphic in all S7. 
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A deep impact into the theory of quadrature domains was the discovery 
by M. Putinar [37] in the mid 1990's of an alternative characterization in 
terms of hyponormal operators. Recall that J. Pincus proved [34] that with 
any bounded linear operator T : H — > H in a Hilbert space H for which the 
self-commutator is positive (i.e., T is hyponormal) and has rank one, say 

[T*,T] = T*T - TT* 0^£eH, 

one can associate a unitary invariant, the so-called principal function. This 
is a measurable function g : C — ► [0, 1], supported on the spectrum of T, 
such that for any z,w in the resolvent set of T there holds 

detfCT.TrX- 1 ) = exp[-L / ^4A% (3) 

2?n Jc (C -*)(C- w ) 

where T u = T—uI. The right hand side in (3) is referred to as the exponential 
transform of the function g. In case g is the characteristic function of a 
bounded set Q we have the exponential transform of 

E n (z,w) = exp[-L / -A- A JL]. (4) 

A central result in Putinar's theory is the following criterion: a domain f2 
is a quadrature domain if and only if the exponential transform of Q is a 
rational function of the form 

En(z,w)= 2 { \»/ = . M>l»l, (5) 

where P and Q are polynomials. In this case Q is the same as the polynomial 
in (2). 

In the present paper we shall unify the above pictures by interpreting 
the exponential transform of a quadrature domain in terms of resultants 
of meromorphic functions on the Schottky double of the domain. To this 
end we need to extend the classical concept of resultant of two polynomials 
to a notion of resultant for meromorphic functions on a compact Riemann 
surface. The introduction of such a meromorphic resultant and the demon- 
stration of its usefulness in several contexts is the main overall purpose of 
this paper. 

The definition of the resultant is natural and simple: given two mero- 
morphic functions / and g on a compact Riemann surface M we define their 
meromorphic resultant as 



where (/) = J2 a i ~ J2 h = — Z -1 ! 00 ) is the divisor of /. This 

resultant actually depends only on the divisors of / and g. It follows from 
Weil's reciprocity law that the resultant is symmetric: 

K(f,g)=K(g,f). 

For the genus zero case the meromorphic resultant is just a cross-ratio 
product of four polynomial resultants, whereas for higher genus surfaces it 
can be expressed as a cross-ratio product of values of theta functions. In 
the other direction, the classical resultant of two polynomials (which can 
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be viewed as meromorphic functions with a marked pole) may be recovered 
from the meromorphic one by specifying a local symbol at the infinity (see 
Section 9). 

It is advantageous in many contexts to amplify the resultant to an elim- 
ination function. With / and g as above this is defined as 

£ f,g( z > w ) = n (f ~ z '9 ~ w), 

where z, w are free complex parameters. Thus defined, £f, g (z, w) is a rational 
function in z and w having the elimination property 

£/, 9 (/(C),<?(0) = (C€M). 

In particular, this gives an explicit formula for the algebraic dependence be- 
tween two meromorphic functions on a compact Riemann surface. Treating 
the variables z and w in the definition of 8f g (z, w) as spectral parameters in 
the elimination problem, the above identity resembles the Cayley-Hamilton 
theorem for the characteristic polynomial in linear algebra. This analogy 
becomes more clear by passing to the so-called differential resultant in con- 
nection with the spectral curves for two commutating ODE's, see for example 
[35]. 

The above aspects of the resultant and the elimination function charac- 
terize them essentially from an algebraic side. In the paper we shall however 
much emphasize the analytic point of view by relating the resultant to ob- 
jects such as the exponential transform (4) and the Fredholm determinant. 
One of the key results is an integral representation of the resultant (The- 
orem 2), somewhat similar to (4). From this we deduce one of the main 
results of the paper: the exponential transform of a quadrature domain Q 
coincides with a natural elimination function on the Schottky double $7 of 

n-. 

Eq(z,w) = 6fj*(z,w). (6) 

Here (/, /*) is a canonical pair of meromorphic functions on S7: / equals 
the identity function on 17, which extends to a meromorphic function on 
the double ft by means of the Schwarz function, and /* is the conjugate of 
the reflection of / with respect to the involution on $7. In Section 8 we use 
formula (6) to construct explicit examples of classical quadrature domains. 

In Section 6 we discuss the meromorphic resultant 7Z(f, g) as a function 
of the quotient 

m = gf m 

Clearly, / and g are not uniquely determined by h in this representation, but 
given h it is easy to see that there are, up to constant factors, only finitely 
many pairs (/, g) with non-zero resultant 1Z(f,g) for which (7) holds. Thus, 
the natural problem of characterizing the total range a{h) of these values 
lZ(f,g) arises. 

Another case of interest is that the divisors of / and g are confined to 
lie in prescribed disjoint sets. This makes 7Z( f, g) uniquely determined by h 
and connects the subject to classical work of E. Bezout and L. Kronecker on 
representations of the classical resultant TZ po \(f,g) by Toeplitz-structured 
determinants with entries equal to Laurent coefficients of the quotient h(z). 
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The 60's and 70's brought renewed interest to this area in connection with 
asymptotic behavior of truncated Toeplitz determinants for rational gen- 
erating functions (cf. [3], [12], [15]). This problem naturally occurs in 
statistical mechanics in the study of the spinspin correlations for the two- 
dimensional Ising model (see, e.g., [6]) and in quantum many body systems 



One of the general results for rational symbols is an exact formula given 
by M. Day [12] in 1975. Suppose that h is a rational function with simple 
zeros which is regular on the unit circle and does not vanish at the origin 
and infinity: ordo h < 0, ordoo h < 0. Then for any N > 1: 



where p, r^, Hi are suitable rational expressions in the divisor of h. An 
accurate analysis of these expressions reveals the following interpretation of 
the above identity in terms of resultants: 



where the (finite) sum is taken over all pairs (/, g) satisfying (7) such that 
g is normalized by g(oo) = 1 and the divisor of zeros of g coincides with the 
restriction of the polar divisor of h to the unit disk: (g) + = (h)- fl D. 

In the above notation, the equality (8) can be thought of as an identity 
between the elements of cr(h) with a prescribed partitioning of the divisor. 
In Section 6.1 we consider resultant identities in the genus zero case in 
general, and show that there is a family of linear relations on a(h). These 
identities may be formally interpreted as a limiting case (for N = 0) of 
the above Day formula (8). Moreover, our resultant identities are similar 
to those given recently by A. Lascoux and P. Pragacz [32] for Sylvester's 
double sums. On the other hand, by specializing the divisor h we obtain a 
family of trigonometric identities generalizing known trigonometric addition 
theorems. Some of these identities were obtained recently by F. Calogero 
in [7, 8]. For non-zero genus surfaces the situation with describing a(h) 
becomes much more complicated. We consider some examples for a complex 
torus, which indicates a general tight connection between resultant identities 
and addition theorems for theta-functions. 

Returning to (6) and comparing this identity with determinantal repre- 
sentation (3) we find it reasonable to conjecture that one can associate to 
any compact Riemann surface an appropriate functional calculus for which 
the elimination function becomes a Fredholm determinant. In Section 7 
we demonstrate such a model for the zero genus case. We show that the 
meromorphic resultant of two rational functions is given by a determinant 
of a multiplicative commutator of two Toeplitz operators on an appropriate 
Hardy space. There are interesting similarities between our determinan- 
tal representation (cf. formula (56) below) of the meromorphic resultant 
and the tau-function for solutions of some integrable hierarchies (see, for 
instance, [43]). 

Further aspects of the meromorphic resultant discussed in the paper 
are interpretations in terms of potential theory, in Section 5, and various 
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cohomological points of view, e.g., an expression of the resultant in terms of 
the Serre duality pairing (subsections 6.3 and 6.4). In Section 4 we give an 
independent proof of the symmetry of the resultant using the formalism of 
currents, and also derive several integral representations. Section 3 contains 
the main definitions and other preliminary material, and in Section 2 we 
review the polynomial resultant. 

The authors are grateful to Mihai Putinar, Emma Previato and Yurii 
Neretin for many helpful comments and to the Swedish Research Council 
and the Swedish Royal Academy of Sciences for financial support. This re- 
search is a part of the European Science Foundation Networking Programme 
"Harmonic and Complex Analsyis and Applications HCAA" . 

2. The polynomial resultant 

The resultant of two polynomials, / and g, in one complex variable 
is a polynomial function in the coefficients of /, g having the elimination 
property that it vanishes if and only if / and g have a common zero [54]. 
The resultant is a classical concept which goes back to the work of L. Euler, 
E. Bezout, J. Sylvester and A. Cayley. Traditionally, it plays an important 
role in algorithmic algebraic geometry as an effective tool for elimination of 
variables in polynomial equations. The renaissance of the classical theory of 
elimination in the last decade owes much to recent progress in toric geometry, 
complexity theory and the theory of univariate and multivariate residues of 
rational forms (see, for instance, [19], [49], [52], [10]). 

We begin with some basic definitions and facts. In terms of the zeros of 
polynomials 

mm n n 

f(z) = /mJJ(z - <k) = ^2fiz\ g(z) =9nYl(z-Cj) = ^2gjZ J , (9) 

i=l i=0 j=l j=0 

the resultant is given by the Poisson product formula [19] 

m n 

n pol (f,9) = /Zll( fl rc;) = &f[g{ai) = (-1)™^ II /to)- (io) 

i,j i=l 3=1 

It follows immediately from this definition that 7Z po \(f,g) is skew-symmetric 
and multiplicative: 

K pol (f,g) = (-i) mn ^ P oi(ff,/), K P oi(/i/2,<7) = Kpoiih, g)Kp6L(f 2, g). 

(11) 

Alternatively, the resultant is uniquely (up to a normalization) defined 
as the irreducible integral polynomial in the coefficients of / and g which 
vanishes if and only if / and g have a common zero. 

All known explicit representations of the polynomial resultant appear 
as certain determinants in the coefficients of the polynomials. Below we 
briefly comment on the most important determinantal representations. The 
interested reader may consult the recent monograph [19] and the surveys 
[10], [49], where further information on the subject can be found. 

With /, g as above, let us define an operator S : V n © V m — ► V m +n by 
the rule: 

S(X,Y) = fX + gY, 
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where Vk denotes the space of polynomials of degree < k — 1 (dim'Pfc = k). 
Then 



^poi(/,5) = det 



(h 



fn 



9o 



fo 9n 



\ 



90 



(12) 



V fm 9n/ 

where the latter is the Sylvester matrix representing S with respect to the 
monomial basis. 

An alternative method to describe the resultant is the so-called Bezout- 
C ay ley formula. For deg / = deg g = n it reads 

TZp i(f,g) = det(/3 i j)o<i i j<„-i, 

where 

f(z)g(w) - f{w)g(z) 



w 



= E 

i,j=0 



Z W 



(13) 



is the Bezoutian of / and g. The general case, say deg / < deg g, is obtained 
from (11) and (13) by completing f(z) to z k g(z), k = degg — deg/. 

Other remarkable representations of the resultant are given as deter- 
minants of Toeplitz-structured matrices with entries equal to Laurent co- 
efficients of the quotient h(z) = ^7jy- These formulas were known al- 
ready to E. Bezout and were rediscovered and essentially developed later 
by J. Sylvester and L. Kronecker in connection to finding of the greatest 
common divisor of two polynomials (see Chapter 12 in [19] and [4]). 

Recently, a similar formula in terms of contour integrals of the quotient 
h(z) has been given by R. Hartwig [28] (see also M. Fisher and R. Hartwig 
[15]). In its simplest form this formula reads as follows. With / and g as 
in (9), we assume go = g(0) / 0. Then for any N > n, the polynomial 
resultant, up to a constant factor, is the truncated Toeplitz determinant for 
the symbol h(z): 

ft P oi(/,<?) = /r N 9 Tl+N dett m , N (h), (14) 
hkz k is the Taylor development of the quotient around 



where h(z 
z = and 



E 



oo 

k=0 • 



( 



l m,N 



(h) 



h 

""m 
h m +i 



h 



-ra- 
th,. 



K 
hm— 



N+l 
N+2 



\ 



hm+N- 



) 



\ Kn+N-l 

and = for negative k. 

The determinant det t m ^{h) is a commonly used object in theory of 
Toeplitz operators. For instance, the celebrated Szego limit theorem (see, 
e.g., [6]) states that, under some natural assumptions, det^o,iv(^) behaves 
like a geometric progression. Exact formulations will be given in Section 7.1, 
where the above identity is generalized to the meromorphic case. 

It is worth mentioning here another powerful and rather unexpected ap- 
plication of dett mi Ar(/i), the so-called Thorn- Porteous formula in the theory 
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of determinantal varieties [18], [20, p. 415]. We briefly describe this identity 
in the classical setup. Consider an n x m (n < m) matrix A with entries a^j 
being homogeneous forms in the variables x±, . . . ,x k of degree pi + qj (for 
some integers pi, qj). Denote by V r the locus of points in F k at which the 
rank of A is at most r. Then, thinking of pi, qj as formal variables, one has 

degV r = dett m _ rjn _ r (c), > c k z = n -. 

k=0 lli=iU Pi*) 

We mention here also a differential analog of the polynomial resultant 
in algebraic theory of commuting (linear) ordinary differential operators. A 
key observation goes back to J.L. Burchnall and TW. Chaundy and states 
that commuting ordinary differential operators satisfy an equation for a 
certain algebraic curve, the so-called spectral curve of the corresponding 
operators (see [36] for a detailed discussion and historical remarks). The 
defining equation of the curve is equivalent to the vanishing of a determinant 
of a Sylvester-type matrix. This phenomenon was a main ingredient of the 
modern fundamental algebro-geometric approach initiated by I. Krichever 
[31] in the theory of integrable equations. By using the Burchnall-Chaundy- 
Krichever correspondence between meromorphic functions on a suitable Rie- 
mann surface and differential operators, E. Previato in [35] succeeded to get 
a pure algebraic version of the proof of Weil's reciprocity. 

All the determinantal formulas given above fit into a general scheme: 
given a pair of polynomials one can associate an operator S in a suitable 
coefficient model space such that 1Z po \{f,g) = detS. On the other hand, 
none of the models behaves well under multiplication of polynomials. This 
makes it difficult to translate identities like (11) into matrix language. One 
way to get around this difficulty is to observe that (14) is a special case 
of the Szego strong limit theorem for rational symbols [15] and to consider 
infinite dimensional determinantal (Fredholm) models instead. We sketch 
such a model in Section 7 below. 

3. The meromorphic resultant 

3.1. Preliminary remarks. For rational functions with neither zeros 
nor poles at infinity, say 

m n 

i=i j=i j 

(\,Hj^0 and all Oj, b{, Cj,dj distinct) it is natural to define the resultant as 
In other words, 

m n , , m n 

w. 9) = n n • iFt = n ik«> ^ m 

. , ■ , u i u i u i °i . . . . 
i=i j=i j ■> i=ij=i 

where (a, b, c, d) := " fzf is the classical cross ratio of four points. 
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Note that (nonconstant) polynomials do not fit into this picture since 
they always have a pole at infinity, but the polynomial resultant can still be 
recovered by a localization procedure (see Section 9). Notice also that the 
above resultant for rational functions actually has better properties than the 
polynomial resultant, e.g., it is symmetric (lZ(f,g) = TZ(g,f)), homogenous 
of degree zero and it only depends on the divisors of / and g. The resultant 
for meromorphic functions on a compact Riemann surface will be modeled 
on the above definition (16) and contain it as a special case. 

3.2. Divisors and their actions. We start with a brief discussion 
of divisors. A divisor on a Riemann surface M is a finite formal linear 
combination of points on M, i.e., an expression of the form 

m 

D = J2mau (18) 
i=i 

a.i G M, rii E Z. Thus a divisor is the same thing as a 0-chain, which acts 
on 0-forms, i.e., functions, by integration. Namely, the divisor (18) acts on 
functions if by 

„ m 

( D ,<P} = / ^ = Vn^(ai). (19) 
Jd i=i 

From another (dual) point of view divisors can be looked upon as maps 
M — > Z with support at a finite number of points, namely the maps which 
evaluate the coefficients in expressions like (18). If D is a divisor as in (18) 
we also write D : M — > Z for the corresponding evaluation map. Then 
D = J2aeM D(a)a. The degree of D is 

m 

deg D = Y,ni = D ( a )- 

i=l aeM 

and its support is 

suppD = {a € M : D(a) ^ 0}. 

If / : M — > P is a nonconstant meromorphic function and a £ ¥ then 
the inverse image f~ 1 (a), with multiplicities counted, can be considered as 
a (positive) divisor in a natural way. The divisor of / then is 

(/) = rHo) - rV). (20) 

If / is constant, not or oo, then (/) = (the zero element in the Abelian 
group of divisors). 

Recall that any divisor of the form (20) is called a principal divisor. In 
the dual picture the same divisor acts on points as follows: 

(/)(«)= ord a (/), 

where ord a (/) is the integer m such that, in terms of a local coordinate z, 

f{z)=c m {z-a) m + c m+1 {z - a) m+l + ... with c m + 0. 

By ord/ we denote the order of /, that is the cardinality of / _1 (0). 

Divisors act on functions by (19). We can also let functions act on 
divisors. In this case we shall, by convention, let the action be multiplicative 
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rather than additive: if h = h(u±, . . . , Uk) is a function and D±, . . . , D k are 
divisors, we set 

h(D l ,...,D k )= 11 h(a u ...,a k ) D ^- D ^\ (21) 

ai,...,a k eM 

whenever this is well-defined. Observe that this definition is consistent with 
the standard evaluation of a function at a point. Indeed, any point a £ M 
may be regarded simultaneously as a divisor D a = a. Then h(a±, . . . , a p ) = 
h(D ai , . . . , D ap ). In what follows we make no distinction between D a and a. 

With branches of the logarithm chosen arbitrarily (21) can also be writ- 
ten 

h(Di,.. .,D p ) = exp (Di ® . . . ® Dp, log /i). 

When Di, i = 1, . . . ,p are principal divisors, say Di = (gi) for some mero- 
morphic functions gi, the definition (21) yields 

h(( gi ), (g p )) = [] h(a u . . . ,a p) ord ai ( ffl )-ord ap ( ffp )_ 

ai,...,a p gM 

3.3. Main definitions. Let now /, g be meromorphic functions (not 
identically and oo) on an arbitrary compact Riemann surface M and let 
their divisors be 

(/) = r\o) - rv) = zr=i * - EL 6 - 

(g)=g- 1 (0)-g- 1 (^) = ^ =i c J -^ =i d J . 

At first we assume that (/) and (g) are "generic" in the sense of hav- 
ing disjoint supports. In view of the suggested resultant (16) for rational 
functions the following definition is natural. 

Definition 1. The (meromorphic) resultant of two generic meromorphic 
functions / and g as above is 

W, g ) = 9 ((/)) = fl = = exp((/),log 9 >. (23) 

In the last expression, an arbitrary branch of logg can be chosen at each 
point of (/). 

Elementary properties of the resultant are multiplicativity in each vari- 
able: 

K(f 1 h,g)=K(f 1 ,g)K(f 2 ,g), K(f,g 1 g 2 )=K(f,g 1 )Tl(f,g 2 ). 

An important observation is homogeneity of degree zero 

K(af,bg)=K(f,g) (24) 

for a, b G C* := C \ {0}. The latter implies that lZ(f, g) depends merely on 
the divisors (/) and (g). 

Less elementary, but still true, is the symmetry: 

K(f,g)=K(g,f), (25) 
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i.e., in the terms of the divisors 

TTfiK) = TT f( c j) 

This is a consequence of Weil's reciprocity law [55], [20, p. 242]. In Sec- 
tion 4 we shall find some integral formulas for the resultant and also give an 
independent proof of (25). 

If, in (21), some of the divisors are principal then the resulting action 
h may be written as a composition of the corresponding resultants. For 
instance, for a function h of two variables we have 

h((f), (g)) = K u {f{u),K v {g{v), h(u, «))), (26) 

where 1Z U denotes the resultant in the u-variable. 

Remark 1. The definition of meromorphic resultant naturally extends 
to more general objects than meromorphic functions. Indeed, of / we need 
only its divisor and g may be a fairly arbitrary function. We shall still 
use (23) as a definition in such extended contexts. However, there is no 
symmetry relation like (25) in general. See e.g. Lemma 4. 

When, as above, (/) and (g) have disjoint supports TZ(f, g) is a nonzero 
complex number. It is important to extend the definition of TZ(f, g) to 
certain cases when (/) and (g) do have common points. 

Definition 2. A pair of two meromorphic functions / and g is said to 
be admissible on a set A C M if the function a — > ord a (</) ord a (/) is sign 
semi-definite on A (i.e., is either > on all A or < on all A). If A = M 
we shall simply say that / and g is an admissible pair. 

It is easily seen that the product in (23) is well-defined as a complex 
number or oo whenever / and g form an admissible pair. 

Clearly, any pair of two meromorphic functions whose divisors have no 
common points is admissible (we call such pairs generic). Another important 
example is the family of all polynomials, regarded as meromorphic functions 
on the Riemann sphere P. It is easily seen that any pair of polynomials is 
admissible with respect to an arbitrary subset A C P. 

The following elimination property is an immediate corollary of the def- 
initions. 

Proposition 1. Let two nonconstant meromorphic functions f, g form 
an admissible pair on M. Then 1Z(f,g) = if and only if f and g have a 
common zero or a common pole. In particular, lZ(f, g) = if f and g are 
polynomials. 

3.4. Elimination function. We have seen above that the meromor- 
phic resultant of two individual functions is not always well-defined (namely, 
if the two functions do not form an admissible pair). However one may still 
get useful information by embedding the functions in families depending on 
parameters, for example by taking the resultant of / — z and g — w. We shall 
see in Section 8.3 that such resolved versions of the resultant have additional 
analytic advantages. 
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Let z, w £ C be free variables. The expression 

£(z, w) = £f, g (z, w) = U(f - z,g -w), 

if denned, will be called the elimination function of / and g. 

Theorem 1. Let f and g be nonconstant meromorphic functions without 
common poles. Then the elimination function is well defined everywhere 
except for finitely many pairs (z,w), and it is a rational function of the 
form 

£{z,w) = , 
P(z)R(w) 

where Q, P, R are polynomials, and 

p( Z )= n r h= n ( w -9(b)). 

Proof. Note that a linear transformation f f — z keeps the polar 
locus unchanged. Thus the elimination function lZ(f — z, g — w) is well- 
defined for all pairs (z,w) such that f~ 1 (z)ng^ 1 (oo) = g~ 1 (w)nf~ 1 (oo) =0. 
Let (z, w) be any such pair. Then applying the symmetry relation (25) we 
obtain 

(g-w)(f-Hz)) _ (/-*)(<?-») 
1 ' ' (g-w)(f-i(w)) (/-z)( & -i(oo))- 
Let /, g have orders m and n, respectively, as in (22), and let {fr 1 } de- 
note the branches of f~ x . Then spelling out the meaning we find, using that 
the symmetric functions of {g(ff 1 (z))} are single- valued from the Riemann 
sphere into itself, hence are rational functions, that 

m 

(g-w)(f-\z)) = U(9(f-\z))-w) = (-ir(w m +R 1 (z)w m - 1 +---+R m (z)), 

i=i 

where the Ri(z) are rational. Similarly, 

(g - w)(f-\oo)) = (-l) m (w m + r lW m ~ l + ■■■ + r m ), 

where the r^ are constants. 

With the same kind of arguments for (f—z)(g~ 1 (w)) and (/— z)(g~ 1 (oo)) 
we obtain 

£(z , _ w m + R^w™- 1 + ■■■ + R m {z) _ z n + + ■ • • + P n {w) 
^ Z,W) ~ w m + r 1 w m - 1 -\ Vr m ~ z n +p 1 z n - 1 -\ V p n 

Clearing the denominators (in the numerators) yields the required state- 
ment. □ □ 

Important, and useful in applications, is the following elimination prop- 
erty of the function £f j9 (z,w). Let us choose ( € M arbitrarily and insert 
z = /(C), w = g(() into £f, g (z, w). Since the functions f — z and g — w then 
have a common zero (namely at Q this gives, by Proposition 1, that 

£/, 9 (/(C),<7(O)=0 (C€M). 

In particular, 

Q(f,g) = o, 

i.e., we have recovered the classical polynomial relation between two func- 
tions on a compact Riemann surface (see [14], [16], for example). 
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3.5. Extended elimination function. We have seen that the elimi- 
nation function is well-defined for any pair of meromorphic functions without 
common poles. One step further, linear fractional transformations allow us 
to refine the definition of elimination function in such a way that it becomes 
well-defined for all pairs of meromorphic functions. 

Namely, let / and g be two arbitrary meromorphic functions and consider 
the function of four complex variables: 

£(z, w;z ,wo) = £f t g(z,w;zo,w ) = 111 -j — -> - — — )■ ( 27 ) 

Let us choose arbitrary the pair (z,Zo). Then we have for divisor: 
( fZz ) = f~ l ( z ) ~ f l { z o)- It is easy to see that the resultant in (27) 
is well defined for any quadruple (z, w; Zq,Wq) with 

[g-\w)Ug-\w )} n \j-\z) U f-\z )] = 0. (28) 

The set X of all (z, w; zq,wq) such that (28) holds is a dense open subset of 
in C 4 . 

Applying then an argument similar to that in Theorem 1, we find that 
the right hand side in (27) is a rational function for (z, w\ zq, wq) € X. We 
call this function the extended elimination function of / and g. 

We have the cross-ratio-like symmetries £(z, w; zq,wo) = £(zo,w ; z, w), 

and 

£(z,w ;z ,w) = — -. 

£{z,w;zo,w ) 

In the case when the elimination function £^ g (z,w) is well-defined we 
have the following reduction: 

£(z,w)£(z ,w ) Q(z,w)Q(z ,w ) 



£(z,w;z ,w ) 



£(z,w )£(z ,w) Q{z,wq)Q{zq,w)' 



with Q as in Theorem 1. 

In the other direction, the ordinary elimination function, if well-defined, 
can be viewed as a limiting case of the extended version. Indeed, it follows 
from null-homogeneity of the meromorphic resultant that 

/ / — z g — w \ 

£(z, w -z , W0 ) = n^ T - WQ , T -^-j, 



and therefore that 



lim £(z, w; zq, u>o) = £(z, w). 



There are still cases when the elimination function is not defined or is 
trivial while its extended version contains information. To illustrate this, 
let us consider a meromorphic function / of order n and let g = f. Then a 
straightforward computation reveals that 

ci \ ( z-zq w- w \ n 

£ f f {z,w;z ,w Q ) = = {z,w,z ,w ) , 

11 \Z — Wq W — Zq ) 

where (z, w, zq, wq) is the cross ratio. 
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3.6. The meromorphic resultant on surfaces with small genera. 

On the Riemann sphere P the resultant reduces to a product of cross ratios 
(17) and the symmetry relation (25) becomes trivial. Note that the cross 
ratio itself may be regarded as the meromorphic resultant of two linear 
fractional functions. 

Prom a computational point of view, evaluation of the meromorphic 
resultant on P is similar to the evaluation of polynomial resultants. Indeed, 
for any admissible rational functions given by the ratio of polynomials, / = 
/1//2 and g = 51/52, one finds that 

n( f ,g) = / ( oor d ~ ( ^(oor d ~ (/) • |4(HIHtH ( 29 ) 

l<n?6i{h,92) ><-pol{j2,gi) 

The latter formula combined with formulas in Section 2 expresses the mero- 
morphic resultant in terms of the coefficients of the representing polynomials 
of / and g. For example, since each resultant in (29) is a Sylvester determi- 
nant (12), 

Ttpoii fi, 9 j) = detS(fi,gj) = det S ij: 
the resulting product amounts to 

H(f,g) = /(oo) ord -^5(oo) ord -^ • det(S^ 1 S 11 S^ 1 S 22 ). 

In Section 7 we give another, more invariant, approach to the represen- 
tation of meromorphic resultants via determinants (see also Section 7.2 for 
the exponential representations of 7Z(f,g)). 

Now we spell out the definition of the resultant in case of Riemann 
surfaces of genus one. Consider the complex torus M = C/L T , where L T = 
Z + rZ is the lattice formed by r G C, Imr > 0. A meromorphic function 
on M is represented as an L T -periodic function on C. Let 



0(C) = MC)= J2 e" i(fc2r+fe(1+ 



r+20) 

k=—oo 



be the Jacobi theta-function. Then any meromorphic function / on M is 
given by a ratio of translated theta-functions: 



and a necessary and sufficient condition that such a ratio really defines a 
meromorphic function is that the divisor is principal, i.e., by Abel's theorem, 
that 

m 

Y^{a t - h) G L. (30) 
i=i 

With / as above and g similarly with Cj and dj, ^j=i( c j — 4/') e * ne 
following representation for the meromorphic resultant on the torus holds: 



w,«7)=nn 



6(cj - ai)9(dj - bj) 
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4. Integral representations 

4.1. Integral formulas. We shall derive some integral representations 
for the meromorphic resultant, and in passing also give a proof of the sym- 
metry (25), Weil's reciprocity law. Let /, g be nonconstant meromorphic 
functions on a compact Riemann surface M of genus p > and recall (23) 
that the resultant can be written 

ft(/,<7)=exp((/),log<7). 

We assume that the divisors (/) and (g) have disjoint supports. Since (/) 
is integer- valued and different branches of log g differ by integer multiples of 
27ri it does not matter which branch of log g is chosen at each point of (/). 
However, our present aim is to treat logg as a global object on M, in order 
to interpret ((/), log g) as a current acting on a function and to write it as 
an integral over M. 

First of all, to any divisor D can be naturally associated a 2-form current 
fiD (a 2-form with distribution coefficients), which represents D in the sense 
that 



( D ,<P) = f= ^A/i D 
Jd Jm 

for smooth functions <p. With D = Y^ n iO<i this fj,p is of course just 

fJ-D = $Ddx Ady = s ^n i 5 ai dx Ady, (31) 

where 5 a is the Dirac delta at the point a and with respect to a local variable 
z = x + iy chosen (only 5 a dxAdy has an invariant meaning). When D = (/) 
we have the following formula. 

Lemma 1. If f is a meromorphic function, then fi^ = ^ d( y ) in the 
sense of currents. 

PROOF. In a neighbourhood of a point a with ord a (/) = m, i.e., 

f(z) = c m (z - a) m + c m+1 (z - a) m+1 + . . . , c m + 0, 

in terms of a local coordinate, we have y = (j^ + h{z))dz with h holomor- 
phic. Hence, 

fdf\ d ( m \ 

d\ — ] = -p— ( h hiz) dz Adz = mir5 a dz Adz = 2^\rab a dx A dy, 

V / J dz \z - a J 

from which the lemma follows. □ □ 

Next we shall make log / and log g single- valued on M by making "cuts" . 
Let o:i,. . . , a p , Pi,. . ■ , (3 P be a canonical homology basis for M such that 
each /3k intersects once from the right to the left (k = 1, . . . ,p) and no 
other crossings occur. We may choose these curves so that they do not meet 
the divisors (/) and (g). 

Since the divisors (/) and (g) have degree zero we can write 

(/) = &Yf, id) = dig 

where 7/, j g are 1-chains. We may arrange these curves so that there are 
no intersections and so that they are contained in M \ (a\ U • • • U (3 P ). 
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Now, it is possible to select single- valued branches of log / and log g in 

M' = M \ (7/ U 7 S U ai U • • • U (3 P ). 

Fix such branches and denote them Log/, Log g. Then Log / and Log g are 
functions, defined almost everywhere on M, and Log g is smooth in a neigh- 
bourhood of the support of (/) and vice versa. In particular, ((/),Logg) 
and {(g), Log/) make sense. 

Now using Lemma 1 and partial integration (with exterior derivatives 
taken in the sense of currents) we get 

K(f,g) = exp((/),Log5) = exp[ / // (/) ALogfif] 

Jm 

r 1 f r,df\ t n r 1 f df IT 

= ex P[7r^/ ALog^r] = exp[— / — AdLogg]. 

27T1 J M / 27T1 J M / 

In summary: 

Theorem 2. Zei / and g be two meromorphic functions on a compact 
Riemann surface whose divisors have disjoint supports. Then 

K{f,g) =exp[— / —AdLogg]. 
2vri J M f 

In particular, for generic z, w, 

1 f df 

£f,g( z > w ) = ex P^ j M J^~z AdL °s(9-w)}. 

It should be noted that the only contributions to the integrals above 
come from the jumps of Log g (and Log (g—w) respectively), because outside 
this set of discontinuities the integrand contains dz A dz = as a factor. 

4.2. Symmetry of the resultant. We proceed to study dLog in de- 
tail. Let first a, b be two points in the complex plane and 7 a curve from b 
to a such that dj = a — b (formal difference). Then, with a single- valued 
branch of the logarithm chosen in C \ 7, 

z — a dz dz z — a 

CtLog — = ; r l[d Arg Thump contribution from 7 

z — b z — a z — b z — b 

dz dz „ . ,„ , N 
2TTidH 7 (z). 



z — a z — b 

Here dH 1 is the 1-form current supported by 7 and defined as the (distri- 
butional) differential of the function which in a neighbourhood of any 
interior point of 7 equals +1 to the right of 7 and zero to the left. Thus 
dH^ is locally exact away from the end points. The function Hj cannot be 
defined in any full neighbourhood of a or b. On the other hand, dH-y is taken 
to have no distributional contributions at a and b. One easily checks that 
this gives a current which represents 7 in the sense that 



dH-y A r 



/7 J M 

for all smooth 1-forms r. Taking r of the form dip gives 



/ 

J A 



d(dHj) A ip = / dHj A dip = dip = tp. 

M JM Jf Jd-y 
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Thus the 0-chain, or divisor, $7 is represented by d(dH 1 ). We can write this 
also as d(dHy) = /ig 7 , where no is defined in (31). Note in particular that 
dH^ is not closed, despite the notation. 

If 7 and a are two curves (1-chains) which cross each other at a point c, 
then it is easy to check (and well-known) that 

dHj A dH a = ±<5 C dx Ady, 

with the plus sign if a crosses 7 from the right (of 7) to the left, the minus 
sign in the opposite case. For the curves a±, . . . , fa in the canonical homology 
basis, the forms dH ai , . . . , dHp p are closed, since the curves are themselves 
closed. 

Now we extend the above analysis to Log / in place of Log |5f • 111 
addition to the jump across 7/ (an arbitrary 1-chain in M \ (a± U . . . U fa) 
with djf = (/)) we need to take into account possible jumps across the 
ak, (3k- In order to reach the right hand side of from the left hand side 
within M' one just follows fa. The increase of Log / along this curve is 
fp y 7 , hence this is also the jump of Log / across from the left to the 
right. With a similar analysis for the jump across fa one arrives at the 
following expression for dLogf: 

This means that 7/ needs to be modified to the 1-chain 

v 

a f=lf + ^2( w ' md P k (f) ■ a k - wind afc (/) • fa), 
k=i 

where, for a closed curve a in general, wind Q (/) stands for the winding 
number 

Notice that daj = djf = (/) and that now Log / can be taken to be single- 
valued analytic in M \ suppoy. The above can be we can summarized as 
follows. 

Lemma 2. Given any meromorphic function f in M there exists a 1- 
chain <jf having the property that do$ = (/), log/ has a single-valued 
branch, Log/, in M \supp<7/ and the exterior differential of Log f, regarded 
as a 0- current in M with jumps taken into account, is 

dLog / = - 2nidH a . 
J 

Since y = the lemma combined with Theorem 2 gives the following 
alternative formula for the resultant. 

Corollary 1. With notations as above 

K(f,g)=exp(- ! ^AdH ag )=exp[ f. (32) 

J M J Jag J 
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In the corollary cr/ may be replaced by any 1-chain 7 with d^f = (g), 
because this will make a difference in the integral only by an integer multiple 
of 2vri. 

Next we compute 

dhogf A dhogg = {% - 2mdH a .) A (— - 2mdH a ) 
J 1 9 

= % AdLogg + dLogf A — + {2mfdH af A dH a 
J 9 

The integral of d Log f Ad Log g = <i(Log f Ad Log g) over M is zero because 
M is closed, and the integral of the last member, (2n\) 2 dH (Jf A dH ag , is 
an integer multiple of (27ri) 2 . Therefore, after integration and taking the 
exponential we get 

exp[r^r / %AdLogg + ^- [ dLog/A^] = l. 

This proves the symmetry: 

Corollary 2. Let f and g be two meromorphic functions on a closed 
Riemann surface with disjoint divisors. Then 

K(f,g)=K(g,f). 

Remark 2. This symmetry is also a consequence of Weil's reciprocity 
law [55] (see Section 9 for further details), and may alternatively be proved, 
in a more classical fashion, by evaluating the integral in Cauchy's formula 
f dM , Log / A dLogg = (cf. [20, p. 242]). It is also obtained by directly 
evaluating the last integral in (32). 

Remark 3. If the divisors of / and g are not disjoint but /, g still form 
an admissible pair, then both lZ(f,g) and 7Z(g,f) are either or oo, hence 
the symmetry remains valid although in a degenerate way In this case, and 
more generally for nonadmissible pairs, Weil's reciprocity law in the form 
(77) (in Section 9) contains more information. 

By conjugating g one gets the following formula for the modulus of the 
resultant in terms of a Dirichlet integral. 

Theorem 3. Let f and g be two meromorphic functions on a compact 
Riemann surface whose divisors have disjoint supports. Then 

|2 „_„„r 1 f d f » d 9i 



\K(f,g)\' = exp[— / f Af ]. (33) 
^ri J m J 9 

Proof. By Lemma 2 we have 

1 1 df ct~ df d~ 

—dLogfAdLogg = — 4 + 4 Adfl*. -dH a A-^- -2mdH a AdH a 
2vri 27ri f g J 9 

Integrating over M and taking the exponential yields, in view of (32), the 
required formula. □ □ 
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5. Potential theoretic interpretations 

5.1. The mutual energy and the resultant. We recall some poten- 
tial theoretic concepts (see, e.g., [41] for more details). The potential of a 
signed measure ("charge distribution") jj, with compact support in C is 

U»(z) = -J log|*-C|dM(0- 

The mutual energy between two such measures, fx and v, is (when defined) 

v) = - J J log \z - C| dfi(z)dv(Q = J U» dv = Ju u d», 

and the energy of \x itself is I(n) = I(fJ,,fi). In case J dv = J d/j, = the 
above mutual energy can after partial integration be written as a Dirichlet 
integral: 

7(/i, v) = — [ dU»A *dU I/ , (34) 
2vr J 

where * is the Hodge star. 

If K C C is a compact set then either I(jx) = +oo for all (i > with 
supp/i C K, J dfx = 1, or there is a unique such measure for which I(ji) 
has a finite minimum value. In the latter case fi is called the equilibrium 
distribution for K because its potential is constant on K (except possibly 
for a small exceptional set), say 

= 7 (const) on K. 

The logarithmic capacity of K is defined as 

cap(K) = e~ 7 = e~ IM . 

(If I(fj) = +oo for all \i as above then cap (K) = 0). 

Now let us think of signed measures as (special cases of) 2-form cur- 
rents. Then, for example, (31) associates to each divisor D in C the charge 
distribution \i = fiD- In particular, for any rational function / of the form 
f( z ) = 111=1 fEf 1 we have the charge distribution 

m m 

v = Hf) = 6a i dx A dy ~ ^ dx A dy > 

1=1 1=1 

the potential of which is = — log |/|. 

One point we wish to make is that the resultant of two rational functions, 
/ and g, relates in the same way to the mutual energy. In fact, with \i = fi(f) 
and v = ^{g), 

\K(f,g)\ 2 = exp[((/),log 5 ) + {(f), T^g)] = e 2 ^^ = e ~ 2 /^ = e ~ 2I ^\ 
hence 

I(H,v) = -log\K(f,g)\. (35) 

The Dirichlet integral (34) for I(/j,, v) essentially gives the link between (35) 
and (33). 
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5.2. Discriminant. Recall that the (polynomial) discriminant Dis po i(/) 
is a polynomial in the coefficients of / which vanishes whenever / has a mul- 
tiple root. In case of a monic polynomial f(z) = YYiLi( z ~ a i) we nave 

Dis po i(/) = (-l) 11 ^ n pol (/,/') = H(a t - aj) 2 . 

i<j 

Thus the discriminant is the square of the Van der Monde determinant. 

The discriminant can be related to a renormalized self-energy of the 
measure \i = The self-energy itself is actually infinite because point 

charges always have infinite energy. Formally: 

„ m 

I(jjl) = J U»d f i= ((/),- log |/|) = - log J] \ ai - aj \ (=+oo). 

M=l 

The renormalized energy 1(h) is obtained by simply subtracting off the 
infinities I(S ai ), i.e., the diagonal terms above: 

T([i) = -logjjjai -dj\ = -logJJ|aj -cij\ 2 = -log|Dis po i(/)|. 

Thus, |Dis P oi(/)| = e" 1 ^) . Here J dfx = deg/ = m, and after normaliza- 
tion (there are m(m — 1) factors in Dis po i(/)) it is known that the transfinite 
diameter 

i 

dooiK) = lim max |Dis po i(/)| , 

ra^oo deg/=m 

equals the capacity: doo(K) = cap (K). 

Notice also that the discriminant may be regarded as a renormalized 
self-resultant 7£ po i (/,/): 

ft P oi(/, /) = Y[(ai ~ aj) IC ^ n Dis P oi(/) = JJ(°* ~ ( 36 ) 

We can use the same renormalization method to arrive at a definition of 
discriminant in the rational case. Let / be a rational function 

f< a = IM = n£i(*-°t) 
n ) f2(z)-ig =1 (z-b i y 

Then applying the scheme in (36) gives 
\f (ai - bj)(bi - aj) 

mai-aJllQH-bj) (37) 

r enorm „ . , . s _ _ /Vpol(Jl, JjJ ><-po\(j2, J 2) 

1S( ): ~ "^ P oi(/i,/ 2 )^ pol (/2,/i)' 

The corresponding renormalized energy of fx = is 



I (jj) = - log 
which yields 



Ui,j(ai-bj)Yli,j(bi-aj) 
iDis(/)| = e- ? M. 



log|Dis(/)| 



20 



BJORN GUSTAFSSON AND VLADIMIR G. TKACHEV 



We note that the definition (37) of Dis(/) is consistent with the so-called 
characteristic property of the polynomial discriminant [19, p. 405]. Namely, 
one can easily verify that the meromorphic resultant of two rational functions 
can be obtained as the polarization of the discriminant in (37), that is 

Dis(/<?) 



Dis(/)Dis( 5 ) 



5.3. Riemann surface case. Much of the above can be repeated for 
an arbitrary compact Riemann surface M. For any signed measure /ionM 
with J M dfx = there is potential U^, uniquely defined up to an additive 
constant, such that 

-d * dU^ = 2irfi. 

Here /x is considered as a 2-form current (/i may actually be an arbitrary 2- 
form current with (/x, 1) = 0, and then will be a 0-current; the existence 
and uniqueness of follows from ordinary Hodge theory, see e.g. [20, 
p. 92]). 

The mutual energy between two measures as above can still be defined 

as 



/( /U , v ) = j V* dv = Ju u 



d\x 



and (34) remains true. Similarly, (35) remains valid for fi = /«(/), v = [i( g y 
Thus 

\n(f,g)\ = e- I ^\ 

It is interesting to notice that this gives a way of defining the modulus 
of the resultant of any two divisors of degree zero: if deg D\ = deg D 2 = 
with supp D\ n supp D2 = then one naturally sets 

\K(D 1 ,D 2 )\ = e- 1 ^ !'^. 

It is not clear whether there is any natural definition of TZ(D±, D2) itself, 
except in genus zero where we have (17). Directly from the definition (23) 
we can however define lZ(D,g) = g(D) for D a divisor of degree zero and g 
a meromorphic function. 

6. The resultant as a function of the quotient 

6.1. Resultant identities. In previous sections we have considered 
the resultant as a function of two meromorphic functions, / and g, say. 
Sometimes, however, it is possible and convenient to think of the resultant 
as a function of just one function, namely the quotient h = jj. In general, 
part of the information about / and g is lost in h, hence some additional 
information has to be provided. 

For instance, if / and g are two monic polynomials, then formula (14) 
in its simplest form, when N = n, reads 

ftpoi(/,5) = det t m , a (h). 

Another example is if the divisors of / and g are confined to lie in 
prescribed disjoint sets: given any set U C M then among pairs /, g with 
supp(/) C U, supp(g) C M \ U, the resultant TZ{f,g) only depends on A 
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Integral representations for 1Z(f,g) in terms of only f/g and U will in such 
cases be elaborated in Section 6.2 (Theorem 4). 

In the remaining part of this section we shall pursue a further point of 
view. Suppose that the divisors of / and g are not necessarily disjoint but 
that / and g still form an admissible pair. In general we have, with h = f/g, 

ord h < ord / + ord g, 

and it is easy to see that lZ(f, g) = if and only if this inequality is strict (be- 
cause strict inequality means that at least one common zero or one common 
pole of /, g cancels out in the quotient f/g)- 

Now start with h and consider admissible pairs f,g with h = f/g and 
such that 

ord/i = ord/ + ordc/. (38) 

In general there are many such pairs /, g and by the above 1Z(f, g) ^ 
for all of them. The question we want to consider is whether there are 
any restrictions on which values TZ( f, g) can take. At least in the rational 
case there turns out be such restrictions and this is what we call resultant 
identities. 

Let d > 1 and 

*w = n^f. < 39 > 

i=l 

Let C™ denote the set of all increasing length-m sequences (h, . . . ,i m ), 
1 < i\ < . . . < i m < d. For two given elements I,Je C™ define 

Then all the solutions /, g of (38), up to a constant factor (which by (24) is 
inessential for the resultant), are parameterized by 

f(z) = h u (z), g(z) = klJ ^ = 1 

h(z) h rjl (z) 

where the prime denotes complement, e.g., I' = {1, . . . , d} \ I. 

The main observation of this section is that the resultants TZ(f, g) satisfy 
a system of linear identities. An extended version of the material below with 
applications to rational and trigonometric identities will appear in [27]. 

Proposition 2. LetO <m <d and J <G C™. Then 

£ n(h u , i/hrj.) = n ( h Ji> yh Jip ) = i. (40) 

Proof. We briefly describe the idea of the proof. Denote by A and B 
the two Van der Monde matrices with entries {a{~ 1 ) and (&!" 1 ), 1 < i, J < d, 
respectively. Let / = {h, ■ ■ ■ ,i m } an d J = {ji, ■ ■ ■ , j m }- Then one can 
readily show that 

H{hu, l/h P j,) = (-If det Ajj det(A" 1 ) /J , (41) 

where n = YlT=i(^s + 3s)- Here A = AB -1 and Ajj (resp. (A~ x )/j) denotes 
the minor of A (resp. A -1 ) formed by intersection of the rows i £ I and 
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the columns j G J. Hence the required identities follow from (41) and the 
Laplace expansion theorem for determinants. □ □ 

In the simplest case, d = 2, m = 1, (40) amounts to the characteristic 
property of the cross-ratio: 

(a, b, c, d) + (a, c, b, d) = 1. 

The resultants in (40) appear also in the so-called Day's formula [12] for 
the determinants of truncated Toeplitz operators. Let h be a function given 
by (39) such that |6j| ^ 1 for all i, and let J = {j : \bj\ > 1}. 

Introduce the Toeplitz matrix of order N 

( h h-i ... h x _ N \ 
h\ ho ... h 2 -N 



t N (h) 



(42) 



y /iJV-l ■ ■ ■ h J 

where h k = ± J e- ike h{e i8 )d9 are the Fourier coefficients of h on the unit 
circle. Then, in our notation, Day's formula reads 

dett N {h) = nhu,l/h r j,) ■ h$j,(0), (43) 

where m denotes the cardinality of J and N > 1. Notice that formal sub- 
stitution of N = with to(h) = 1 into (43) gives exactly the statement of 
Proposition 2. 

Remark 4. Taking double sums in (40) (over all I,J E Cf) we get 
quantities which occur also when computing subresultants (see, e.g., [32]). 
Recall that the (scalar) subresultant of degree k is the determinant of the 
matrix obtained from the Sylvester matrix (12) by deleting the last 2k rows 
and the last k columns with coefficients of /, and the last k columns with 
coefficients of g. In a different context, the subresultants are determinants of 
certain submatrices of the Sylvester matrix (12) which occur as successive 
remainders in finding the greatest common divisor of two polynomials by 
the Euclid algorithm [50] . 

The identities (40) have beautiful trigonometric interpretations. Take 

■Pr z- e 2iak -i— r z — e 2ic ' 
= 11 z _ e 2ib fc ' 9{z) = 11 J3^id7- 
k=l 1=1 

Then one easily finds that 

hence a direct application of (40) gives the following. 
Corollary 3. Let d > 2 and J € C™. Then 

x - Ujj' s Hai ~ bj>) Ui',j s[n ( b j ~ a i') = (u) 

i n.i,i> sin K - a i') rijj' sin ( b j ~ b j) ~ ' { ' 

where the sum is taken over all subsets I G C™ and the product over i G I, 
i> G /', j G J, f G J'. 
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For example, specializing by taking bj = § + aj in (44) one gets identities 
in the spirit of those given recently in [7] , [8] . 

There are also analogues of Proposition 2 for the complex torus M = 
C/L T . For these one has to take into account the Abel condition (30). 
Although we have not been able to find complete analogues of the rational 
resultant identities, one particular case is worth mentioning here. Notice 
that the minimal possible value of d in order for a meromorphic function 
h( z ) = Y\i=i t° split into two non-constant meromorphic functions, 

i.e. h = f /g, is d = 4. One can readily show that any such function may be 
written as 

_ 4>{z - ZQ,ai)(j>(z - z ,a 2 ) 
4>{z - z ,bi)4>(z - z ,b 2 ) ' 
where 4>((, a) = 9(( — a)9(( + a). We additionally assume that a\ ± d 2 L 
and b\ ± b 2 L. Then all non-constant solutions of (38) are given by 

„. . 6(z,a,j) , , 6(z,bji) 
<p{z,bj) (p{z,ai>) 

where {k,k'} = {1,2}. Hence 



Pij ■■= K(f,g) 



9(ai — bji)0(a,i + bj')9(a,i' — bj)9(di> + b 



1 2 



3) 



9(di — dii)9{di + di>)9(bj — bj')Q(bj + bj')j ' 
and there only two different values of pif. 

Ci := Pn = P22, ^2 := P12 = P21- 
Using the famous addition theorem of WeierstraB 

=9( a - c)9(d + c)9(b - d)9(b + d) - 6 {a - b)9(d + b)9(c - d)9(c + d) 
-9( a - d)9(d + d)9(b - c)9(b + c), 
one finds that (with appropriate choices of signs) 

±y/h±y/& = l, (45) 

or more adequately: (1 — £i) 2 + (1 - ^) 2 = 2^i^2- 

The identity (45) may be generalized to functions of the kind 

h ( z > = 11 IT/ iTT" 

However the problem of description of the range of 7Z( f, g) in (38) for general 
meromorphic functions h on C/L T remains open. 

6.2. Integral representation of TZjj. Let us now turn to the situation 
of having a preassigned set U C M and consider resultants lZ(f,g) for 
meromorphic functions / and g with supp(/) C U, supp(g) C M \ U. It 
is easy to see that for such pairs 7Z(f, g) only depends on the quotient h = 
f/g. Indeed, this is obvious from the fact (see (24)) that the resultant only 
depends on the divisors: under the above assumptions the divisors of / and 
g are clearly determined by h and U. 

To make the above in a slightly more formal we may define 1Z(Di, D 2 ) 
for any two principal divisors D±, D 2 having, e.g., disjoint supports. For 
any divisor D, let Djj denote its restriction to the set U and extended by 
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zero outside U (thus with D = YlaeM D(a)a, Djj = Yla&u D(a)a). Then in 
the situation at hand we can write 

K(f, g) = K((f), (<?)) = K((h)u, (h)u - (h)), 
which only depends on h and U. This motivates the following definition. 

Definition 3. For any set U C M and any meromorphic function h on 
M such that (h)u is a principal divisor we define 

n u {h) = n{(h) u ,(h) u -(h)). 

It is easy to check that 

TZu(h) = TZ M \u(h). 
We shall consider the symmetric situation that 

M = U U T U V, 

where U, V are disjoint nonempty open sets and V = dU = dV. We provide 
r with the orientation of dU. By the above, with / and g meromorphic on 
M, supp(/) C U, supp(<7) C V and h = f /g we have 

Hu{h)=K v {h)=K{f,g). 

Note that the function h is holomorphic and nonzero in a neighbourhood of 
r, h G 0*(T), and that it is uniquely defined by its values on T. Our aim is 
to find an integral representation for lZjj(h) in terms only of the values of h 
on r. 

The problem of decomposing a given h G 0*(T) into functions / G 
0*(V), g G 0*(U) with h = f/g is a special case of the second Cousin 
problem. By taking logarithms we shall reduce it, under symplifying as- 
sumptions, to the corresponding additive problem, which is the first Cousin 
problem. For the latter we have the following simple criterion for solvability. 

Lemma 3. Let M = U U T U V be as above. Necessary and sufficient 
condition for a function H G 0(F) to be decomposable as 

H = H + — H- onT 

with H + eO(U),H-€ 0{V) is that 

J H A to = for all cue O 1 ' (M). 

When the decomposition exists the functions H± are unique up to addition 
of a common constant (more adequately: a function in 0(M)). 

The lemma is well-known and can be deduced for example from the Serre 
duality theorem. We shall just remark that "explicit" representations of H± 
can be given in terms of a suitable Cauchy kernel: 

H±(z) = ^~j T H((Wz, C; z , Co) d( 

the plus sign for z G U, minus for z G V . The kernel &(z, £; zq, Co) is, in the 
variable z, a meromorphic function with a simple pole at z = ( and a pole 
of higher order (depending on the genus) at z = Co- I n the variable C it is 
a meromorphic one-form with simple poles of residues plus and minus one 
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at ( = z and ( = zq respectively; zq and Co are fixed but arbitrary points, 
zq 7^ Co- in the case of the Riemann sphere, <3?(z, (; zq, Co) d( 1S the ordinary 
Cauchy kernel 

*(z,C;«d,Co)dC = ? 3-- 7 3-, (46) 

hence does not involve Co- in the the case of higher genus the point Co is 
really needed. We refer to [40] for the construction of the Cauchy kernel in 
general. 

Theorem 4. Let M = UUTUV with U connected and simply connected, 
and let h be meromorphic on M without poles and zeros on T. Assume in 
addition that 

1 f dh 

-r = (47) 

2vri J T h V ' 

and that 

Log/iAi<J = for allu G O 1 ' (M) (48) 

T 

(the previous condition guarantees that a single-valued branch of log h exists 
on T). Then (h)u is a principal divisor and 

K u (h)=exp [^r^d(Log/i)_ A (Log /»)+]. 

Remark 5. Ideally (48) should be replaced be the weaker condition that 
there exists a closed 1-chain 7 on M such that 

J Log h A u = 2vri J u for all lo G 1 '°(M). (49) 

In fact, this turns out to be exactly, by Abel's theorem, the necessary and 
sufficient condition for (h)jj to be a principal divisor. However, (49) would 
lead to a more complicated formula for IZu(h). Note that (48) is vacuously 
satisfied in the case M = P, which will be our main application. Condition 
(47) says that the divisor (h)u has degree zero. 

Proof. We first prove that (h)u is a principal divisor. Using the nota- 
tion of Lemma 2 we make Log h into a single- valued function on all of M by 
making cuts along a 1-chain ah such that dah = (h). Since Log/i is already 
single- valued on T, ah can be chosen not to intersect T. Thus ah consists 
of two disjoint parts, ahP\U and ah H V. The terms of ah containing the 
curves ai, . . . , f3 p will appear in cr^ n V because U is simply connected. 

Now, for all uj G O 1,0 (M) we have by (48) and Lemma 2 

= — f Log h A uj = — !- / cCLog h A u = — !— / f 2mdH a , ) Aco 

2vri 7 r 2vri 7^ 2vri 7^ V h ) 



1 /" d/i 
2m J v h 



Alo — I dH ah Alo = — dH ahCi u Auj = — I lo. 
u Jm ' Ja h nu 



By Abel's theorem this implies that d(ah D U) = {h)jj is a principal 
divisor (condition (49), in place of (48), would have been enough for this 
conclusion) . 

The divisor (h)u being principal means that (h)u = (f) for some / 
meromorphic on M. Setting g = f/h we have supp(/) C U, supp(g) C V 
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and h = f/g. It follows that IZu(h) = lZ(f,g), hence to prove the theorem 
it is by Theorem 2 enough to prove that 

[ d(Logh)_ A(Log/i)+ = I ^f- Ad Log 5. 
Jr Jm J 

To that end we shall compare two decompositions of dLogh = * on T: 
from Lemma 3 we get 

dLog h = d(Log h)+ — d{Log on V 

with (Log/i)+ G 0(U), (Log/i)_ G 0(V), while h = f/g gives 

dh df dg 



h 



f 



9 



on r, 



where df/f G O 1 ' ^), dg/g G O l >°{U). 
It follows that 



f 



+ d(Log h). 



dg_ 
9 



+ d(Log h). 



on r 



and that the left and right members combine into a global 1-form ujq G 
lfi (M). Thus 



d(Log/i). 



u - — in I/, 



<i(Log = u; 



dg_ 
9 



in [/. 



In the simply connected domain U we may write loq = d<p for some 
(p G 0{U) and also ^ = dLogg (dH CTg = in U because a g can be chosen 
to be ah fl F; similarly oy can be chosen to be <7/j PI U). It follows after 
integration and adjusting 99 by a constant that 

(Log h) + = ip - Log g in U. 

Now we finally obtain 

^d(Log/i)_ A (Log/i)+ = ^(u; - y) A (<£>-Log#) = - jf y A (99 - Log 5) 

= / — AdLog^- / (dLog/i + dLogs-) A 99 = / — AdLogg, 
J Jr Jm J 



as desired. 



□ 



□ 



Remark 6. Under the assumptions of the theorem, the solution of the 
second Cousin problem of finding /, g such that h = f/g on T is given by 



/ = exp 
9 = exp 



df 
f 



I 

fdg 
J 9 _ 



exp 
exp 



{u — d(Log/i)_ 
J (uj — d(Log h). 



in V, 
in U 



(indefinite integrals), where to G O 1,0 (M) is to be chosen such that J(uj — 
<i(Log/i)_) is single- valued in V modulo multiples of 2ni. 
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6.3. Cohomological interpretations of the quotient. Let us give 
some interpretations of the above material in terms of Cech cohomology. 
Given h G 0*(T), let U±, V\ be open neighbourhoods of U and V, respec- 
tively, such that h G 0*{U\ C\V\). Then {U\, V{\ is an open covering of M, 
and relative to this h represents an element [h\ in i7 1 (M, O*). It is well- 
known [21], [16] that [h] = as an element in # X (M, O*) if and only if h 
is a coboundary already with respect to {Ui,V±}, i.e., if and only if there 
exist / G 0*(Vi) and g G 0*(U X ) such that h = f/g in Ui D Vl. If /i is 
meromorphic in M, then so are / and g. 

Similarly, a function H G O(r) represents an element [H] in H l (M,0), 
and [#] = if_and only if there exist F G 0(E7i), G G 0(Vi) (for some 
Ui dTJ,ViD V) such that H = F - G on T. 

The spaces H l (M,0) and H l (M,0*) are related via the long exact 
sequence of cohomology groups which comes from the exponential map on 
the sheaf level: with e(f) = exp[27ri/] we have 

0->Z-»OAO*->1, 

hence 

-» fl°(M, Z) -»• iJ°(M, O) fl°(M, O*) -» ^(M, Z) -» 
-» ^(M^O) ^ H l (M,0*) H 2 (M,Z) -» 0. 
From this we extract the exact sequence 

-» H l {M,0)/H l (M,Z) A H l {M,0*) A H 2 (M,Z) 0. (50) 

Here c is the map which associates to [/i] G if 1 (M, C>*) its characteristic 
class, or Chern class, and it is readily verified that it is given by 

c ([/i]) = wind r h = ^[ J r ~f^ = de s( h )u- 

If c([/i]) = then [/i] is in the range of e. If T is connected then log fa 
is single-valued on T and the preimage of [h] can be represented by H = 

Log h. However, if T is not connected then the preimage of [h] cannot 
always be represented by a function H G 0(T), one needs a finer covering 
of M than {Ui, V±} to represent it. This is a drawback of the method using 
the decomposition M = U UT UV in combination with the exp-log map 
and explains some of our extra assumptions in Theorem 4. 

Assume nevertheless that the preimage of [h] G H l (M, O*) (with c([h]) = 
0) can be represented by H = ^Log/i G 0(T). Then of course [h] = 
if [H] = as an element in H l (M,0), i.e., if J r H A u = for all to G 
1,O (M). However, what exactly is needed for [h] = is by (50) only that 
[H] G H 1 (M,Z,), and this what is expressed in (49). 

Since, for H G 0(T), [H] = as an element in H\M, O) if and only if 
Jr H A = for all to G O 1 ' (M), the pairing 



descends to a bilinear map 

H°(M, O 1 ' ) x i? 1 (M, 0) -► C. 
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This map is in fact the Serre duality pairing ([44], [21]) with respect to the 
covering {U±, V{\. Versions of the Serre duality with respect to more general 
coverings will be discussed in the next section. 

6.4. Resultant via Serre duality. We now return to the general in- 
tegral formula in Theorem 2, and interpret the exponent Jm y A ciLog g 
directly in terms of the Serre duality pairing, which in general also involves 
a line bundle or a divisor. With a divisor D, the pairing looks 

( , ) Serre : H°(M, 0]f) x H\M, 0- D ) - C, 

between meromorphic (l,0)-forms with divisor > — D and (equivalence 
classes of) cocycles of meromorphic functions with divisor > D. 

In our case, given two meromorphic functions / and g, we choose D > 
to be the divisor of poles of y (or any larger divisor), so that y £ 

T(M, O^f). As for the other factor, log g defines an element, which we 
denote by [S log g], of i/ 1 (M, O-d) as follows. First, with j g as in the be- 
ginning of Section 4.1, choose an open cover {Ui} of M consisting of simply 
connected domains Ui satisfying 

(supp D U supp j g ) PI Ui fl Uj = whenever i ^ j 

(in particular supp j g ndUi = for all i). Second, choose for each i a branch, 
(logg)j, of \ogg in Ui\~/ g . Finally, define a cocycle {{S log g)ij}, to represent 
[<51og<7] eH\M,0- D ), by 

(5 log g)ij = (logg)j - (log g)j \nUi n Uj. 

There exist smooth sections over Ui, vanishing on D, such that 

{S log g)ij = i>i - ipj mUiDUj. (51) 

One may for example choose a smooth function p : M — > [0, 1] which vanishes 
in a neighbourhood of supp D U supp 7 9 and equals one on each UiPiUj, i ^ j 
and define ipi = p(logg)i in Ui. In any case, (51) shows that the ipi satisfy 

dipi = dipj in Ui n Uj , 

so that {dipi} defines a global (0, l)-form dip on M. The Serre pairing is 
then defined by 

(^,[6logg]) Seiie = ^ [ ^A8^. 

It is straightforward to check that the result (mod27ri) does not depend 
upon the choices made, and that it (mod2-7ri ) agrees with J M y A dLogg. 

A variant of the above is to consider the product y A [5 log g] directly as 
an element in i7 1 (M, O 1,0 ), because there is a natural multiplication map 

H°(M, 0]f) x H l {M, 0_ D ) -» H\M, O 1 ' ), 

and use the residue map (sum of residues; see [21], [16]) 

res : H X (M, O 1,0 ) -» C. 

Then one verifies that 

res (— A [S log g\) = / — Adhogg (mod27ri). 

J 27ri JM J 
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In summary we have 

Theorem 5. For any two meromorphic functions f and g 

df df 
K(f,g) = exp((y, [51og 5 ])scrrc) = exp(res(y A [Slogg])). 

The above expressions can be viewed as polarized and global versions 
of the torsor, or local symbol, as studied by P. Deligne, see in particular 
Example 2.8 in [13]. 

7. Determinantal formulas 

7.1. Resultant via Szego's strong limit theorem. In this section 
we show that the resultant of two rational functions on P admits several 
equivalent representations, among others as a Cauchy determinant and as a 
determinant of a truncated Toeplitz operator. We start with establishing a 
connection between resultants and Szego's strong limit theorem. 

Let us apply the results of the previous section to the case when 

M = P, U = B, V = P\B, r = T = <9D, 

and h is holomorphic and nonvanishing in a neighbourhood of T with windT h - 
(equivalent to that log fa has a single- valued branch on T in this case). 
Choose an arbitrary branch, Log h, and expand it in a Laurent series 

oo 

Log h(z) = ^s k z k . 

— oo 

Note that sq is determined modulo 27riZ only and that the s k also are the 
Fourier coefficients of Logh(e ld ): 

s k = (Log h) k = i- J 2 J e - ifce Log h(e ie ) M. (52) 

Then using the Cauchy kernel (46) with zq = oo one gets 

oo oo 

(Log h) + {z) = ^s k z k , (Log h)-{z) = ~Y, s -kZ~ k , 

k=0 k=l 

and d(Log h)-(z) = YlkLi ks-k-^rr- This gives the formula 

oo 

Tlo(h) = exp[y~] ks k s-k\. (53) 
k=i 

In particular, we have the following corollary of Theorem 4. 

Corollary 4. Let f and g be two rational functions with supp(/) C D 
and supp(^) CP\D. Then 

j, oo 

K(f, g) = TJd(-) = exp[^ fcs fc a_ fc ], (54) 
9 fc=i 

where Log ^i^y = SfcL-oo s k^ lke is the corresponding Fourier series. 
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The right member in (54) admits a clear interpretation in terms of the 
celebrated Szego strong limit theorem (see [6] and the references therein). 
Indeed, under the assumptions of Corollary 4, 

f(JS\ 00 

m^) = ^ = J2 h ^ ke e L °°(n 

therefore h naturally generates a Toeplitz operator on the Hardy space 
tf 2 (B): 

T(h) : d> -> P+(/^), 

where <fr € H 2 (H>) and P + : L 2 (T) — > i? 2 (D) is the orthogonal projection. 
Denote by i(/i) the corresponding (infinite) Toeplitz matirx 

t(h)ij = hi-j, i,j > 1 

in the orthonormal basis {e lke }k>o- 

Then the Szego strong limit theorem says that, after an appropriate 
normalization, the determinants of truncated Toeplitz matrices dettjv(/i) 
(defined by (42)) approach a nonzero limit provided h is sufficiently smooth, 
has no zeros on T and the winding number vanishes: windx(/i) = (see [6], 
[47]). 

To be more specific, under the assumptions made, the operator T(l/h)T(h) 
is of determinant class (see for the definition [47, p. 49]) and 



lim e -N(Logh)o dett N {h) = expV fe(Log/i) fc (Log/i)_ fe = detT(l/h)T(h), 

(55) 



Af->oo 

fc=l 



where (Log/i)^ = Sk are defined by (52). Thus IZn(h) = det T (I /h)T(h). 

We have the following determinantal characterization of the resultant 
(cf. (3)). 

Proposition 3. Under assumptions of Corollary 4, the multiplicative 
commutator 

r( 5 )r(/)- 1 r( 5 )- 1 r(/) 

is of determinant class and 

K(f,g) = det nL )T (l) = det[T(f)~ 1 T(g)T(f)T(g)~ 1 ] 

= lim (4t-K) -dettjv(-) C5fil 

oo 

= exp^/c(Log /i) fe (Log/i)_ fc . 
fc=i 

PROOF. In view of Corollary 4, it suffices only to establish that the 
operator determinants and the limit in (56) are equal. Assume that / and 
g are given by (15). Then 

W 9 M 9(0) M 1 - % /i i - i 
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Expanding the logarithm 

, . . f(oo) v^ T 1 — di/z v^ T 1 — z/d 



in the Laurent series on unit circle \z\ = 1 we obtain: (Log/i)o = Log 
and 

! f E^iK" fc -6r*). if fc <° 

(Log/i) fc = 



*' ! >:; if /•• > 

By the assumptions on the zeros and poles of / and g, this yields that 
J2kez I^H( L °g h )k\ 2 < oo. By the Widom theorem [56] (see also [47, p. 336]) 
we conclude that T(h)~ 1 T(h) — I is of trace class. Therefore the Szego 
theorem becomes applicable for h(z). Inserting the found value (Log/i)o 
into (55) we obtain 

lim (4r^) ■debt N (h)=detT(l/h)T(h). 
N ^°° V/i 00 )/ 

It remains only to show that 

T(i/h)T(h) = T(/)- 1 r( 5 )T(/)r( 5 )- 1 . 

In order to prove this, notice that by our assumptions g,l/g £ H 2 (0) with 
sup zeD \g(z)\ < oo, and f(l/z) e H 2 (B) with inf 2eD \f(l/z)\ > 0. Thus 
h(z) = f(z)/g(z) is the Wiener-Hopf factorization (see, for example, [47], 
Corollary 6.2.3), therefore T(h) = T(f)T(l/g) = T(f)T(g)~ 1 . Similarly we 
get T(l/h) = T(/)^ 1 T( 5 ) and desired identity follows. □ □ 

7.2. Cauchy identity. A related expression for the resultant for two 
rational functions is given in terms of classical Schur polynomials. Namely, 
the well-known Cauchy identity [48, p. 299, p. 323] reads as follows: 



ran 



= ^S\(a)S\{c) = ex.p^2kp k (a)p k (c). (57) 
i=ij=i- aiC i a k=i 

Here A = (Ai, A2, . . . , A&, . . .) denotes a partition, that is a sequence of non- 
negative numbers in decreasing order Ai > A2 > . . . with a finite sum, 

, , / \j+m—j\ , , 1 \j+m— j\ 

S X ( X ) . s x ( Xl , X2 , ...) = = det % - ^ 

det(x])i<ij< m 11 (xi-Xj) 

l<i<j<m 

stands for the Schur symmetric polynomials and 

I m 1 n 

i=i j=i 

are the so-called power sum symmetric functions. 

Note that the series in (57) should be understood in the sense of formal 
series or the inverse limit (see [33, p. 18]). But if we suppose that 

\ ai \ < 1, \cj\ < 1, Vt, j, (58) 

then the above identities are valid in the usual sense. 
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Let us assume that (58) holds. In order to interpret (57) in terms of the 
meromorphic resultant, we introduce two rational functions 

m n 

/(z) = II(l-^), g{z) = lkl-zci 

We find 



i=i z j=i 



yV ' 1 = 1 j = l 

and by comparing with (57) we obtain 



K{f,g) = exp[-^ kp k (a)p k (c)]. (59) 



/ 



k=i 

By virtue of assumption (58), supp(/) G D and supp(^) GP\D, which 
is consistent with Corollary 4. One can easily see that (59) is a particular 
case of (54). 

8. Application to the exponential transform of quadrature 

domains 

8.1. Quadrature domains and the exponential transform. Abounded 
domain 17 in the complex plane is called a (classical) quadrature domain [1], 
[42], [46], [26] or, in a different terminology, an algebraic domain [53], if 
there exist finitely many points Zi G U and coefficients Cj G C (i = 1, . . . , N, 
say) such that 

N 

h dxdy = 2. Cih(zi) (60) 
n i=i 

for every integrable analytic function h in £1. (Repeated points Zj are allowed 
and should be interpreted as the occurrence of corresponding derivatives of 
h in the right member.) 

An equivalent characterization is due to Aharonov and Shapiro [1] and 
(under simplifying assumptions) Davis [11]: £2 is a quadrature domain if 
and only if there exists a meromorphic function S(z) in Q (the poles are 
located at the quadrature nodes zi) such that 

S(z) =z for z G dn. (61) 

Thus S(z) is the Schwarz function of dft [11], [46], which in the above case 
is meromorphic in all of £1. 

Now let be an arbitrary bounded open set in the complex plane. The 
moments of Q are the complex numbers: 

z m z n dxdy. 

Recoding this sequence (on the level of formal series) into a new sequence 
b mn by the rule 

oo , oo 

EOmn -i / V""^ "mil x 1111^1 

____ = i_ eX p(_ ^ ____), |z|,m»i, 

m,n=0 m,n=0 
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reveals an established notion of exponential transform [9], [38], [23]. More 
precisely, this is the function of two complex variables defined by 



E n (z,w) = ex Pi^[J^ 



(- z (-W 

It is in principle defined in all C 2 , but we shall discuss it only in (C \ Q) 2 , 
where it is analytic/antianalytic. 

For large enough z and w we have 

b 

Eq(z,w) = 1— ^ _mn_ 

m,n=0 

Remark 7. The exponential transform admits the following operator 
theoretic interpretation, due to J.D. Pincus [34]. Let T : H — > H be a 
bounded linear operator in a Hilbert space H, with one rank self-commutator 
given by 

[T*,T] = T*T - TT* = £<g>£, 

where £ G H, £ ^ 0. Then there is a measurable function g : C — > [0, 1] with 
compact support such that 

detlTXrr 1 ^- 1 ] = exp[-^ / /°^ AdC " ], (62) 

where T u = T — ul. The function </ is called the principal function of T. 
Conversely, for any given function g with values in [0, 1] there is an operator 
T with one rank self-commutator such that (62) holds. 

Let O be an arbitrary bounded domain. In [37] M. Putinar proved that 
the following conditions are equivalent: 

a) Q is a quadrature domain; 

b) £1 is determined by some finite sequence (a mn )o< m ,n<A r ; 

c) for some positive integer TV there holds 

det(6 mra )o<m,n<Af = 0; 

d) the function Eq(z,w) is rational for z,w large, of the kind 

T-l / X Q(Z,W) 

Eq(z,w) = -=== , (63) 
P{z)P(w) 

where P and Q are polynomials; 

e) there is a bounded linear operator T acting on a Hilbert space H, 
with spectrum equal to f2, with rank one self commutator [T*,T] = 
£ <8> £ (£ £ H) and such that the linear span V/fc>o T* k £ is finite 
dimensional. 

When these conditions hold then the minimum possible number N in 
b) and c), the degree of P in d), and the dimension of Vfc>o-^* fc £ m e ) an 
coincide with the order of the quadrature domain, i.e., the number N in 
(60). For Q, see more precisely below. 

Note that Eq is Hermitian symmetric: Eq(w,z) = Eq(z,w) and multi- 
plicative: if fli and ^2 are disjoint then 

(z, w) = E ni (z, w)E n , 2 (z, w). (64) 
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As \w\ — > oo one has 

En(z,w) = l-^K n (z) + 0(-^) (65) 
w \w\ z 

with z G C fixed, where Kq(z) = / stands for the Cauchy trans- 

n 

form of fL On the diagonal w = z we have Eq(z,z) > for z € C \ $7 
and 

lim z) = 

for almost all zo € 30 (see [23] for details). Thus the information of dfl is 
explicitly encoded in Eq. 

It is also worth to mention that 1 — Eq(z,w) is positive definite as a 
kernel, which implies that when f2 is a quadrature domain of order N then 
Q(z,w) admits the following representation [24]: 

N-l 

Q(z, w) = P(z)P(w) - Pk(z)P^), 

k=0 

where deg Pk = k. 

In the simplest case, when Q = D(0,r), the disk centered at the origin 
and of radius r, the Cauchy transform and the Schwarz function coincide 

2 

and are equal to and 

E o(o,r)(z,w) = 1 -. (66) 

v ' ' zw 

8.2. The elimination function on a Schottky double. Let Q be a 

finitely connected plane domain with analytic boundary or, more generally, 
a bordered Riemann surface and let 

m = n = nudnun 

be the Schottky double of fi, i.e., the compact Riemann surface obtained 
by completing f2 with a backside with the opposite conformal structure, the 
two surfaces glued together along d£l (see [14], for example). On Q there 
is a natural anticonformal involution : f2 ^ f2 exchanging corresponding 
points on Q and O and having dQ as fixed points. 

Let / and g be two meromorphic functions on Q. Then 

r = (w), <?* = (w)- 

are also meromorphic on O. 

Theorem 6. With 0, 0, /, g as above, assume in addition that f has 
no poles inUU dQ and that g has no poles in QU dfl. Then, for large z, w, 

£f,g{z,w) = exp — / A= — = . 

2m J n f-z g*-w 

In particular, 

n , s r 1 f d f d l 1 

^ m Jn J — z J — w 
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PROOF. For the divisors of / — z and g — w we have, if z,w are large 
enough, supp(/ — z) C 17, supp(g — w) C 17. Moreover, log(g — w) has a 
single- valued branch in 17 (because the image g(17) is contained in some disk 
O(0, R), hence (g — u>)(17) is contained in D(— w, R), hence log(g-w) can be 
chosen single- valued in 17 if \w\ > R). Using that g = g* on 917 we therefore 
get 

1 f df 1 /* dj~ 

£f, g (z,w) =exp[— / AdLog(# - tD)] = exp[— / AdLog(g 

Jn J - z 2m J n f - z 

r 1 f df , . 1 f df 

= ex P -:r^/ 7 ALog(p - = exp -— / A Log 

2vri Janf-z 2vri J dn f - z 

^ f df dg~* . 
= exp[ — : / — ^— A 



2niJ n f-z g* -w 
as claimed. □ □ 



8.3. The exponential transform as the meromorphic resultant. 

Let S(z) be the Schwarz function of a quadrature domain f2. Then the 
relation (61) can be interpreted as saying that the pair of functions S(z) 
and z on Q combines into a meromorphic function on the Schottky double 
Q, = U d£l U of Q, namely the function g which equals S(z) on 17, z on 

n. 

The function f = g* = g o <fr is then represented by the opposite pair: 
z on Q, S(z) on 17. It is known [22] that / and g = f* generate the 
field of meromorphic functions on 17, and we call this pair the canonical 
representation of 17 in Q 

From Theorem 6 we immediately get 

Theorem 7. For any quadrature domain 17 

E n (z,w) = S f j*(z,w) (\z\, \w\ » 1), 

where f, f* is the canonical representation of 17 in 17. 

Here we used Theorem 6 with /(C) = ( on 17, i.e., /|n = id. A slightly 
more flexible way of formulating the same result is to let / be defined on 
an independent surface W, so that / : W — > 17 is a conformal map. Then 17 
is a quadrature domain if and only if / extends to a meromorphic function 
of the Schottky double W (this is an easy consequence of (61); cf. [22]). 
When this is the case the exponential transform of 17 is 

Eq(z,w) = £fj*(z,w), 
with the elimination function in the right member now taken in W. 

Remark 8. If 17 is simply connected one may take W = D, so that 
W = P with involution <p : ( h-> Then / : D —> 17 is a rational function 
when (and only when) 17 is a quadrature domain, hence we conclude that 
Eq(z,w) in this case is the elimination function for two rational functions, 

/(C) and /*(C) = /(V0- to Pi c w ^ De pursued in Section 8.5. 
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In analogy with (27) one can also introduce an extended version of the 
exponential transform: 

tp < \ r 1 f ( \ ( d( d( 

E n (z, w; z ,w ) := exp — / A -= r - -= 

2inJ n \(-z C-z J \(-w C-w , 

One advantage with this extended exponential transform is that it is 
defined for a wider class of domains, for example, for the entire complex 
plane. If the standard exponential transform is well-defined then 

p / x E n (z,w)E n (z ,w ) 

E n (z,w;z ,w ) = — — — — -. 

Eq(z,w )Eq(z ,w) 

In other direction, the standard exponential transform can be obtained from 
the extended version by passing to the limit: 

Eq(z,w) = lim Eq(z,w;z ,wo). 

zo,wo— >oo 

Arguing as in the proof of Theorem 7 we obtain the following general- 
ization. 

Corollary 5. Let Q is a quadrature domain with canonical represen- 
tation f and /*. Then 

E n (z, w; z ,w ) = Sfj* 0, w; z ,w ), 

where £fj*(z,w;zo,wo) is the extended elimination function (27). 

8.4. Rational maps. Now we study how the exponential transform of 
an arbitrary domain in M = P behaves under rational maps. For simplicity, 
we only deal with bounded domains, but this restriction is not essential. It 
can be easily removed by passing to the extended version of the exponential 
transform. 

For domains in general, the exponential transform need not be rational. 
However we still have the limit relation (65). This makes it possible to 
continue Eq at infinity by 

Eq(z, oo) = Eq(oo,w) = £h(oo,oo) = 1. 

Theorem 8. Let Q i} i = 1,2, be two bounded open sets in the complex 
plane and F be a p-valent proper rational function which maps fii onto ^2- 
Then for all z,w € C \ 0,2 

E p 2 {z,w) = E X {{F - z),(F - w)) =K u {F{u) - z,TZ v (F(v) - w^E^v))), 

(67) 

where Ek = En k . (See (21) for the notation.) 
Proof. We have 



(F(C)-z)(F(0-w) 

1^2 *'l 



Let D u denote the divisor of F(() — u. Then 

F '<« -f,o g ,F K) -,) = Eft La) 



no-* <*c ^ r c-a 



THE RESULTANT ON COMPACT RIEMANN SURFACES 



37 



where the latter sum is finite. Conjugating both sides in this identity for 
z = w we get 



F'(() _ DJJ3)_ 
2^ 7_7? ' 



therefore, 



(F(Q - z)(F(C) -w) (C - a)(C - 0) 

By assumptions, F(Q — u is different from and oo for any choice of 
u G C \ 1^2 and ( G fii. Hence suppZ) u C C \ VL\. Thus successively taking 
the integral over £l\ and the exponential gives 

E p 2 (z,w) = J] ^ 1 (a,/3) D ^ a ) D -(' 3 )=i5; 1 ( J D 2 , J D ?i ,), 
which is the first equality in (67). Applying (26) we get the second equality. 



Since the exponential transform is a hermitian symmetric function of its 
arguments, a certain care is needed when using formula (67). The lemma 
below shows that the meromorphic resultant is merely Hermitian symmetric 
when one argument is anti-holomorphic. Indeed, suppose, for example, that 
/ is holomorphic and g is anti-holomorphic, that is g(z) = h(z), where h is 
a holomorphic function. Note that (g) = (h). Therefore 

K(g, f) = f((g)) = f((h)) = h((f)) = gW)) = W7g). 

In summary we have 

Lemma 4. Let f(z) be holomorphic (or anti-holomorphic) and g(z) be 
anti-holomorphic (holomorphic resp.) in z. Then 

n{gJ) = K{U). (68) 

Corollary 6. Under the conditions of Theorem 8, if E\ is rational 
then E\ is also rational. 

Proof. First consider the inner resultant 1Z V (-, ■) in (67). Since E±(u, v) 
and F(v) — w are rational and E\ is hermitian, the resultant is a rational 
function in u and w by virtue of (29) and Sylvester's representation (12) 
(see also Lemma 4). Repeating this for 1Z U (-, •) we get the desired property. 

□ □ 

Remark 9. The fact that rationality of the exponential transform is 
invariant under the action of rational maps is not essentially new. In the 
separable case, that is when Eq 1 is given by a formula like (63), and in 
addition / is a one-to-one mapping, the rationality of Eq 2 was proven by 
M. Putinar (see Theorem 4.1 in [37]). This original proof used existence of 
the principal function (see Remark 7). 
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8.5. Simply connected quadrature domains. Even for quadrature 
domains, Theorem 8 provides a new effective tool for computing the expo- 
nential transform and, thereby, gives explicit information about the complex 
moments, the Schwarz function etc. 

Suppose that is a simply connected bounded domain and F is a uni- 
formizing map from the unit disk D onto £1. P. Davis [11] and D. Aharonov 
and H.S. Shapiro [1] proved that fi is a quadrature domain if and only if F 
is a rational function. The we have (cf. Remark 8). 

Theorem 9. Let F be a univalent rational map of the unit disk onto a 
bounded domain f2. Then 

E n (z,w)=K u (F(u)-z,F*{u)-w) (69) 



where F*(u) = 

PROOF. We have from (66) that En(u,v) = 1 — i. Hence E®(u, •) has 
a zero at i and a pole at the origin, both of order one. Applying (68) we 
find 



TZ v (F(v) - w,E (u,v)) = TZ v (E (u,v),F(v) - w) 



Fjp-w _ F*fr)-w 
F(0) - w F(0) - w ' 



Taking into account the null-homogeneity (24) of resultant and using The- 
orem 8 we obtain (69). □ □ 

Applying (29) can we write the resultant in the right hand side of (69) 
explicitly. 

Corollary 7. Let F(Q = be a univalent rational map of the 

unit disk onto a bounded domain fi, where B is normalized to be a monic 
polynomial. Then 

E n (z, W )=KMB,Bt) n *f 1 ^\ (70) 

T(z)T(w) 

where m = degB, n = max(deg A, deg B) = degF, P t = A — tB, 
T(z) = (F(0) - z) n - m 1l pol (P z ,Bt), 



and -P"(C) = C degP -P(l/C) i s the so-called reciprocal polynomial. 

We finish this section by demonstrating some concrete examples. First 
we apply the above results to polynomial domains. Let, in Corollary 7, 
F(C) = aiC + ■ ■ ■ + a n ( n be a polynomial. Then B = fl» = 1, T(z) = z n and 

P z (() = -z + oiC + • • • + a n C, Pl(C) = a n + ... + ai C~ l ~ wC- 
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This gives the following closed formula. 

Z" 1 



En(z, w) = det 



On 
W 



ai 



On 

z 



-1 2±L 

w 



9± 

z 



-1 



V 



ai 



w 



(71) 



1/ 



A similar determinantal representation is valid also for general rational func- 
tions F. 

For n = 1 and n = 2, (71) becomes 
Eq(z,w) = 1 - xiyi, 
Eq(z,w) = 1 - xiyi - 2x 2 j/2 



2 2 
X 2V2 



where Xj = ai/z and 7/j = aj/-u). 

The determinant in (71), and, more generally, the resultant in (69), has 
the following transparent interpretation in terms of the Schwarz function. 
Suppose that VL = F(H>) for a rational function F and recall the definition 
(61) of the Schwarz function of d£l: S(z) = z, z G dft. After substitution 
z = F(Q, |C| = 1, this yields 



(72) 



S(F(0) = F(C)=F(-) = F*(Q. 

Note that F*(() is a rational function again. Thus the Schwarz function 
may be found by elimination of the variable £ in the following system of 
rational equations: 

(w = F*(C), 
\ z = F(Q, 

where w = S(z). Namely, by Proposition 1 the system (72) holds for some 
( if and only if 

K c (F(Q-z,F*(Q-w) = 0. (73) 

The latter provides an implicit equation for w = S(z) in terms of z. 
Note that the expression on the left hand side in (73) is exactly the expo- 
nential transform Eq(z,w) in (69). In fact, Theorem 7 implies that for any 
quadrature domain Q one has Eq(z, S(z)) = 0. 



9. Meromorphic resultant versus polynomial 

Recall that the meromorphic resultant vanishes identically for polyno- 
mials (considered as meromorphic functions on P) . This makes it natural to 
ask whether there is any reasonable reduction of the meromorphic resultant 
to the polynomial one. Here we shall discuss this question and show how to 
adapt the main definitions to make them sensible in the polynomial case. 

First we recall the concept of local symbol (see, for example, [45], [51]). 
Let f,g be meromorphic functions on an arbitrary Riemann surface M. 
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Notice that for any a 6 M, the limit 

\ord a g 



r a (f,g) : =(-l)-da/ord a51im /(£) 



z^a g(z) wAa f 

exists and it is a nonzero complex number. This number is called the local 
symbol of /, g at a. 

For all but finitely many a we have T a (f, g) = 1. The following properties 
follow from the definition: 

Ta(f,g)r a (g,f) = l, (74) 

multiplicativity 

r a (f,g)r a (f,h)=T a (f,gh), (75) 

and 

Ta(/,<7)° rda/l r a ( 5 ,/ l ) 0rd ^r a (/l,/) 0rd ^ = (_l)°rd a /.ord a9 .ord fc (?6) 

In this notation, Weil's reciprocity law in its full strength states that 
on a compact M, the product of the local symbols of any two meromorphic 
functions / and g equals one: 

X\r a (f,9) = l. (77) 

aeM 

Definition 4. Let a £ M and let / and g be two meromorphic functions 
which are admissible on M \ {a}. Let a = a(() be a local coordinate at a 
normalized such that a(a) = 0. Then the following product is well-defined: 

K*(f,9) = Ta(a, °rV U^T^ 1 (78) 

and is called the reduced (with respect to a) resultant. 
Proposition 4. Under the above assumptions, 

K a (f,9) = (-l) orda/ ord » s • K*(g,f), (79) 

and 

KMif2,9) =K a {h,g)n a {h,g). (80) 
Moreover, if a' is another local coordinate with cr'(a) = 0, then 

KAf,9) = (-n(a',o-)r d »f ^°KM,g)- (81) 

PROOF. Note first TZ c {f,g) vanishes or equals infinity if and only if 
lZ a (g,f) does so. Indeed, let us assume that, for instance, lZ a (f,g) = 0. 
Then it follows from (78) and the fact that r a (-, •) is finite and never vanishes, 
that y(£o)° rde ° (/) = for some £o 7^ a - Hence ordg (/) ordg (g) > 0, and 
/ (Co) ^ ^ = 0. From the admissibility condition we know that the product 
ord^(/) ord^(g) does not change sign on M\{a}, therefore ord^(/) ord^(^) > 
everywhere. Then changing roles of / and g in (78), we get 1Z a {g, f) = 0. 

Thus without loss of generality we may assume that TZ a (f,g) ^ and 
T^-a(f,g) 7^ oo. By virtue of the definition of admissibility we see that the 
product ord^ /ord^ g is semi-definite on M \ {a}, hence 

ord 5 /ord 55 = (£€M\{a}). (82) 
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Since ord a cr = 1, we have by (76) and (74) 

_ f\ord a g 

T a a \^oM a f = r a (9^r daf r a (a,fr d ^ = (-l) ord « ford ^r a (g, f) 



We have 

TMgJ) = T a (f,g)r a (a,fr d ^ fJQ^kW 
K*(f,9) r a (g,f)r a (a,g)o^af n g ^or di (f) 

f ( C\ord^(g) 

- (_-|\ord a /ord a g ( f \ TT J Is; 

= (-ir dafordag r a (f,9)ii(-ir^ foi ^ 9 n(f,g). 

Hence, by virtue of (82) and (77) we obtain 

^-o-(g) /) _ \ord a /ord a g TT ft \ _ t i \ord a /ord a g 

and (79) follows. 

In order to prove (80), it suffices to notice that the right side of (78) is 
multiplicative, by virtue of (75), with respect to /. 

Finally, we notice that by (76): r a (a' , g) Ta {g, a)r a (a, a'f Tda9 = (-l) or( H 
hence 

?M/.g) _ ( Ta(<T',g) \ 0tdaf _ (_ T („> ^ord ofl ord a / 

Ha(f,g) ~W(a,g)) " ( Ta{a > a)) 

and the required formula (81) follows. □ □ 

Now we apply some of the above constructions to the polynomial case. 
On the Riemann sphere, P, we pick the distinguished point a = oo and the 
corresponding local coordinate a(z) = -. Since any two polynomials form 
an admissible pair on C, the corresponding product in (78) is well-defined. 

Let us consider two arbitrary polynomials / and g. Since ord^ /-ord^ g > 
for any point £, we see that 1Z a (f,g) = if and only if / and g have a 
common zero in C. In particular, lZ a (f,g) ^ for coprime polynomials. 

Now let / and g have no common zeros. In the notation of (9) we have 
ordoo g = —n and 

r 00 (a,g) = (-ir lim — = { —^- 
z-oc g(z) g n 

and 

f(z)~ n a m 
Too (f,g) = {-l) nm lim J ) I = (-\) nm ^- 

hence 

m n 

KM, g) = &U g(0 OTddf) = /X II lite - <*) 

£7^00 i=lj=l 

Thus, comparing this with (10), we recover the classical definition of poly- 
nomial resultant. We have therefore proved the following. 
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Corollary 8. Let M = P and a{z) = - be the standard local coordinate 
at oo. Then 

K a (f,g) = K po i(f,g). 

A beautiful interpretation of the product in the right hand side of (78) 
as a determinant is given in a recent paper of J.-L. Brylinski and E. Previato 
[5]. In particular, the authors show that this product is described as the 
determinant det(/, A/gA) of the Koszul double complex for / and g acting 
on A = H (M\{a},O). 
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